Effects of dimensionality and anisotropy on the Holstein polaron 
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We apply weak-coupling perturbation theory and strong- 
coupling perturbation theory to the Holstein molecular crys- 
tal model in order to elucidate the eflFects of anisotropy on 
polaron properties in D dimensions. The ground state en- 
ergy is considered as a primary criterion through which to 
study the effects of anisotropy on the self-trapping transi- 
tion, the self-trapping line associated with this transition, 
and the adiabatic critical point. The effects of dimension- 
ality and anisotropy on electron-phonon correlations and po- 
laronic mass enhancement are studied, with particular atten- 
tion given to the polaron radius and the characteristics of 
quasi-lD and quasi-2D structures. Perturbative results are 
confirmed by selected comparisons with variational calcula- 
tions and quantum Monte Carlo data. 

PACS numbers: 71.38. -fi, 71.15.-m, 71.35.Aa, 72.90. -fy 



I. INTRODUCTION 

Polarons are ubiquitous quasiparticles in deformable 
materials embodying the renormalizing effects of defor- 
mation quanta (phonons) on free carriers. The effects 
that can appear depend on the strength of the electron- 
phonon coupling and on the relative time scales of the free 
electron motion and the relevant host vibrations, this lat- 
ter relationship being subsumed in the common notion 
of adiabaticity. Loosely speaking, at fixed adiabatic- 
ity weakly-coupled polarons may be spread over many 
lattice sites ("large" or "free") while strongly-coupled 
polarons may be highly localized, even essentially com- 
pletely collapsed ("small" or "self-trapped"). Similarly, 
at fixed electron-phonon coupling strength very adiabatic 
polarons may be quite broad, while non-adiabatic po- 
larons may be quite compact. Part of the looseness in 
this characterization has to do with the nature of the 
self-trapping transition, and with dependences on the ef- 
fective dimensionality of the host system. 
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Well-known and widely-invoked results tied to the adi- 
abatic approximation |p]-|l6| suggest that polarons in ID 
should be qualitatively distinct from those found in 2D 
and 3D, and that even the notion of self-trapping should 
take on different meaning in low and high dimensions 
p]-p^. The root of this lies in stability arguments sug- 
gesting that in ID all polaron states should be character- 
ized by finite widths, while in 2D and 3D polaron states 
may have either infinite radii ( "free" states at weak cou- 
pling) or finite radii ( "self-trapped" states at strong cou- 
pling) . The self-trapping transition is thus taken to mean 
the abrupt transition from delocalized "free" states char- 
acterized by the free electron mass to highly-localized 
"self-trapped" states characterized by strongly-enhanced 
effective masses. 

The overall conclusion of this paper, on the other hand, 
consistent with a growing body of independent work 
|l7|-p2| , is that the properties of higher-dimensional po- 
larons are more qualitatively similar in most respects to 
those of ID polarons than they are different, and that 
those distinctions that can meaningfully be drawn are 
only distantly related to the more familiar expectations 
outlined above. 

Central in this subject is the notion of anisotropy, 
which figures particularly strongly in quasi-lD systems 
such as conducting polymers or in quasi-2D systems such 
as high-Tc materials. The nature of self-trapping in 
a quasi-lD system poses particularly potent challenges, 
since depending on one's stance one may reach divergent 
conclusions: the polaron may self-trap, or it may not; it 
may be sharply localized, soliton-like, or free; its mass 
may be the free electron mass, may be weakly renor- 
malized, or may be enhanced by orders of magnitude. 
The resolution of such ambiguities lies not uniquely at 
the interface between ID and 2D systems, but in a gen- 
eral understanding of the role of anisotropy in polaron 
structure, within which the clarification of the nature of 
self-trapping in quasi-lD systems is a byproduct. 

An essential preliminary observation is that the self- 
trapping transition is not, in fact, an abrupt phenomenon 
in any dimension except in the adiabatic limit; at finite 
parameter values the physically-meaningful transition is 
more in keeping with a smooth, if rapid, "crossover" from 
polaron structures characteristic of the weak-coupling 
regime to structures characteristic of the strong-coupling 
regime. The "self-trapping line" describing this transi- 
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tion can be located by criteria sensitive to changes in po- 
laron structure; these may involve physical observables 
such as the polaron ground state energy and effective 
mass, or may rely upon more formal properties less ac- 
cessible to direct physical measurement. Here, we con- 
sider several physical observables at finite parameters as 
well as in asymptotic regimes. Through these, we are 
able to characterize self-trapping in one, two, and three 
dimensions for any degree of anisotropy. 

For the explicit calculations to follow, we use the Hol- 
stein Hamiltonian |^,^ on a D-dimcnsional Euclidean 
lattice 



H — Hkin + Hph + Hint , 
D 



n i=l 



Hki, 

Hph ^fiuj'^ bibn , 

n 

Hint = -ghoj ^ a\ari{b\ 4- 6,-j) , 



(1) 
(2) 
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in which at creates a single electronic excitation in the 

rigid-lattice Wannier state at site r?, and fet creates a 
quantum of vibrational energy in the Einstein oscillator 
at site n. All sums are understood to run over the entire 
infinite, periodic, D-dimensional lattice. Because there 
is no phonon dispersion in this model, and because the 
electron-phonon coupling is strictly local, it is in H^in 
where lattice dimensionality and structure have their 
greatest influence; the Ji are the nearest-neighbor elec- 
tronic transfer integrals along the primitive crystal axes, 
and the ti are unit vectors associated with the primitive 
translations. The above model encompasses all Bravais 
lattices, with the different lattice structures appearing 
only in the relative values of the hopping integrals Ji. 
For simplicity in the following, we use terms appropriate 
to orthorhombic lattices in which conventionally i — x, y, 
or z; however, all results hold for lattices of lower sym- 
metry with appropriate transcription of these labels to 
those of the primitive axes. 

For some purposes in this paper, we qualify and quan- 
tify anisotropy through a vector 



(5) 



whose orientation in a Cartesian system can be used 
to objectively quantify anisotropy. In these terms, an 
isotropic property is one depending only on the modulus 

|J| = (J.J)l/2. 

For other purposes, however, it is convenient to think 
of dimensions being turned "on" or "off" according to 
whether particular Ji are finite or vanishing. Anisotropy 
can then be tuned by varying selected Ji in the inter- 
val (0, J]. In several illustrations to follow, we do this 
sequentially, so that dimensions are "turned on" one by 
one, arriving ultimately at the isotropic i'-dimensional 



case in which Ji = J along all axes. In the following, 
we reserve the unsubscripted scalar symbol J to repre- 
sent common magnitude of all Ji in an isotropic case 
{ J — Ji — \J\/\/D). This manner of tuning dimensional- 
ity does not isolate anisotropy, however, since changing 
one Ji keeping others fixed changes both the orientation 
and modulus of J. 

From either perspective, it is such tuning between di- 
mensions by continuously varying the anisotropy that 
is the physically-meaningful concept in most situations 
characterized as quasi-lD {Jx » Jy, Jz) or quasi-2D 

{Jx, Jy » Jz)- 

Quite apart from such very direct quantifications of 
anisotropy, another quantity that arises naturally in the 
following is the sum of the transfer integrals along each 
axis 



cx Tr Mq ^ 



(6) 



in which Mq is the reciprocal effective mass tensor of 
the free electron (see Eq. ^if.); in the isotropic case, 
J = DJ. 

Using weak-coupling perturbation theory (WCPT) 
identifying the unperturbed Hamiltonian as Hq = Hkin + 
Hph and the perturbation as H' = Hint p^-|27[|, one can 
show that the form of the polaron energy band at weak 
coupling in D dimensions is given by 



E{K)^E\^^^{K) + E\^'^{K) + 0{g''} 



(2) 



where 



cos Kj 



D 



g^n^uj^ dt e-'"^*r[e-2-''*'=°'^'='/o(2J,i), 



(7) 



(8) 



(9) 



in which /„(z) is the modified Bessel function of order n. 

Using strong-coupling perturbation theory (SCPT) fol- 
lowing the Lang-Firsov transformation, identifying the 
unperturbed Hamiltonian as Hq = Hph -\- Hint and the 
perturbation as H' = Hkin) [^,|7-29|, one finds 



E{H) = Efcin) + E^hn) + Ef^in) + O 



J^ 



(10) 



where J is an effective, dressed tunneling parameter that 
may be either comparable to or much smaller than the 
bare J depending on regime, and 



(11) 
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E''g^(K) = —e ^ 2Ji cos Kj 

i=l 

D 



(12) 



D 



— e f{g^)'^^JiJj cos Hi cos Kj 
/(y) = Ei(y) - 7 - ln(7;) 



(13) 
(14) 



where 7 is the Euler's constant and Ei(?;) is the exponen- 
tial integral. 
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FIG. 1. Ground state energies (in units of Tuo) in ID, 
2D, and 3D for J/Tioj = 1. Solid line: Strong-coupling 
perturbation theory through second order. Dashed line: 
Weak-coupling perturbation theory through second order. 
Circles: Quantum Monte Carl data kindly provided by P. 
E. Kornilovitch pl| , p2[ . Triangles: Estimated self-trapping 
points gsT as discussed in Section 
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In Figure |l|, we show results for the ground state en- 
ergy in ID, 2D, and 3D for weak-coupling perturbation 
theory through second order and strong-coupling per- 
turbation theory through second order, together with 
corresponding results of quantum Monte Carlo simula- 
tion 1 21 2^. This comparison shows that in any di- 
mension weak-coupling and strong-coupling perturbation 
theory are both quite good up to a relatively small in- 
terval around the "knee" that is associated with the 
self-trapping transition. The real knee as discernible in 
the quantum Monte Carlo data characteristically falls 
to the strong-coupling side of the intersection (or near- 
intersection) of the WCPT and SCPT curves, but to 
the weak-coupling side of gsr as discussed in Section |^ 
pQ|. (These systematic offsets are symptomatic of the 



smooth nature of the physically-meaningful transition, 
as discussed in the next section.) 

Beyond validating both weak and strong-coupling per- 
turbation theory as used in this paper. Figure |l| also 
holds a message regarding the qualitative character of 
self-trapping in different dimensions. That message is 
that the occurrence of self-trapping is qualitatively simi- 
lar in one, two, and three dimensions, with the primary 
qualitative changes being that the transition systemati- 
cally increases in abruptness and shifts to stronger cou- 
pling as dimensionality is increased. 



II. SELF-TRAPPING PRELIMINARIES 

A central result upon which the following sections build 
is the concept of a self-trapping line as contained in the 
empirical curve 



5st[1] = 1 + y/j/huj 



(15) 



that has been found to accurately characterize the tran- 
sition between the small and large polaron regimes in one 
dimension. This curve was inferred through the applica- 
tion of objective criteria to physical properties such as the 
polaron effective mass ||3l|], ground state energy, kinetic 
energy, phonon energy, electron-phonon interaction en- 
ergy pOl , and electron-phonon correlation function [ p2[ . 
Although this simple construct consistently describes a 
wealth of data drawn from multiple polaron properties 
obtained by our own and independent methods, it is well 
to stress what the above relation is not, since the same 
limitations apply to other constructs we are led to in the 
balance of this paper. 

The self-trapping curves we address do not describe 
a phase transition, nor even the exact location of the 
objectively-determined point of crossover implicit in any 
one physical property. Different physical properties gen- 
erally signal the occurrence of self-trapping at distinct, 
though systematically and tightly clustered points on the 
polaron phase diagram. The empirical self-trapping line 
is not intended to describe any one property exactly, but 
to accurately describe the central trend of clusters of 
transition properties over a large range of the polaron 
parameter space. Self-trapping loci drawn from observa- 
tions of different physical properties thus track the em- 
pirical trend line with their own systematic deviations 
that narrow as the adiabatic limit is approached. 

Clearly, the form of 55t[1] has not been derived from 
first principles. Indeed, one can easily be persuaded from 
approximate descriptions of the problem that the initial 
dependence of f;sT[l] upon J/Hlu is most likely not sin- 
gular, but regular, if perhaps steep ||3^]. Our retention of 
the square root in gsT [1] and its higher-dimensional gen- 
eralizations developed below thus reflects not an assertion 
of singular physical behavior, but merely an economy of 
phenomenology. 
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III. SELF-TRAPPING TRANSITION IN THE 
SCALING REGIME 

The location of the self-trapping line can be estimated 
under the practical assumption that the self-trapping 
transition should lie near the crossover from the weak- 
coupling regime to the strong-coupling regime; here, 
specifically, by the intersection of the ground state energy 
curves as given by the leading orders of each perturbation 
theory 

This is an imperfect assumption, since the errors 
in both weak-coupling perturbation theory and strong- 
coupling perturbation theory increase in absolute terms 
as the transition is approached. Absolute precision in the 
perturbative energy is not required, however, in order to 
accurately locate the transition in parameter space. Han- 
dled carefully, we can expect a WCPT/SCPT crossing 
condition to capture dependences that are asymptotically 
correct in the adiabatic strong-coupling limit, provided 
that relative errors in the appropriate quantities remain 
controlled. Since we are limited to the low orders of per- 
turbation theory, we necessarily depend upon there being 
no unexpected surprises lurking in the higher orders of 
either weak or strong-coupling perturbation theory that 
upset the scaling relationships evident in the leading or- 
ders; that such might, in principle, occur is a caveat, 
however unlikely, that must attach to our arguments. 

To this end, we consider the adiabatic strong-coupling 
limit where the self-trapping line coincides with the adia- 
batic critical point; i.e., where the smooth physical tran- 
sition steepens critically. The composite parameter 



A 



(16) 



appears frequently in discussions of this regime because 
it embodies the essential scaling relationship character- 
izing the adiabatic strong-coupling limit. This dominant 
scaling relationship, ~ J' /hoj, guides our application 
of perturbation theory to the estimation of the location 
of the self-trapping transition. We note that A depends 
not on the modulus of J but on J , the sum of the com- 
ponents Ji, and thus by its very definition A includes a 
dependence on anisotropy. 

The perturbative results ^ - (|l^ can be used to infer 
the expected value of the adiabatic critical point by re- 
taining only those terms that dominate in the adiabatic 
strong-coupling regime; these are terms of comparable, 
leading magnitude when both g and J jhu) are large such 

that g^ ^ J /hoj. The weak-coupling correction -E^^ is 
0{— A} at large J and is thus negligible relative to e'^u, 



which is —2J. The strong-coupling correction E^g^ is 

exponentially small (in g^) relative to E^^q and is thus 
negligible in the adiabatic strong-coupling regime. Of 
the several contributions to E^'g^ appearing in (^3|), only 
the term containing f(2g^) is not exponentially small, 
and of the terms contributing to this non-exponential 



contribution, only a single, dominant term remains non- 
vanishing in the adiabatic strong-coupling regime. Thus, 
combining all non- vanishing terms through second order 
of both weak- and strong-coupling perturbation theory, 
the crossing condition that obtains is 



-2J = -g'^hui - 



\Jl_ 

g^huj 



(17) 



The graphical solution of (HTh is indicated in Figure 




FIG. 2. Ground state energies in the adiabatic strong- 
coupling limit for the isotropic case. Horizontal lines indicate 
the l.h.s. of ([itI ) and curved lines indicate the r.h.s. of (p^. 
Solid: ID. Short-dashed: 2D. Long-dashed: 3D. Symbols in- 
dicate Ac. The chain-dotted curve {—g^) can be viewed as 
corresponding to the D — case. 

In obtaining (|l^, we are using perturbative results 
in extreme limits that may not obviously lie within the 
scope of the retained orders of either perturbation the- 
ory, or in principle may even lie beyond the scope of one 
or the other perturbation theory taken to all orders. Al- 
though the legitimacy of our arguments in this regard is 
beyond the scope of any available proof, it is not unsup- 
ported; that the trends in the true ground state energy 
are consistent with the scaling properties used to obtain 
( p^ is evident in the results of multiple independent non- 
perturbative methods on both sides of the self-trapping 
transition Q. Such studies are necessarily at finite pa- 
rameter values, however, and though confirmatory can- 
not in themselves prove that these trends continue un- 
abated into the adiabatic limit. In Appendix A, we pro- 
vide a discussion of WCPT in particular, showing how it 
is feasible that WCPT may continue to be valid in the 
scaling regime despite what may appear to be essentially 
strong coupling. 

We note that it is the second-order SCPT contribution 
on the r.h.s. of ( pT[ ) that is the crucial element in much 
of the discussion that follows. If one fails to capture this 
contribution to the crossing condition, the self-trapping 
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criterion that results is simply A = 1. This is, in fact, a 
widely-asserted self-trapping condition and is not grossly 
incorrect in many cases; however, considerable structure 
is lost to the casualness with which this estimate is often 
used, and the potential exists for significant quantitative 
errors if applied to the inappropriate regimes. 

In terms of the coraposite parameter A and according 
to the full condition (^_7|), the adiabatic critical point in 
any dimension is given by 



A. 



1± \/l 



(18) 



Of the two roots, it is the larger, root that is the 
physically meaningful one, yielding the isotropic (super- 
scripts "z" ) critical values 



1 r 
2 



l + Vl-D- 



(19) 



Ac[l] = 0.5 , Xl[2] = 0.8536... , A'jS] = 0.9082... , (20) 



Ae[0] = , A^[oo] = 1.0 



(21) 



The (— ) root, besides implying an unmeaningful depen- 
dence of the ground state energy on parameters, would 
imply a Ac that decreases with increasing dimensionality, 
contrary to considerable evidence, including the quantum 
Monte Carlo data shown in Figures |l| and ^. 

The dependence of the adiabatic critical point on 
anisotropy contains interesting structure (see Figure |^): 
In all of ID and in each weakly anisotropic case for D > 1, 
the dependence of Ac (solid line) on the anisotropy is es- 
sentially flat; thus, in the generic case of ordinary bulk 
materials with only modest anisotropies, Ac would ap- 
pear to change significantly with dimensionality but to be 
essentially insensitive to underlying anisotropy. On the 
other hand, in the case of "low-dimensional" materials 
characterized by weak tunneling into one or more trans- 
verse dimensions (e.g., quasi-lD scenarios with one or two 
transverse dimensions, or quasi-2D scenarios with one), 
the weakly-involved dimensions have relatively strong ef- 
fects on Ac: The transition between zero dimensions and 
any higher dimensional case is marked by a jump discon- 
tinuity in the dependence of Ac on any Ji . The transition 
between one dimension and any higher-dimensional case 
is marked by a square-root singularity in dependence of 
Ac on the transverse Ji. The transition between succes- 
sive higher dimensional cases is generically smooth, how- 
ever, with the appearance of any singularity being de- 
pendent on the manner in which dimensionality is tuned 
(see Figure I if.). 
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FIG. 3. Dependence of \st on anisotropy, and related 
curves. Solid line: Ac (adiabatic limit) as anisotropy is 
tuned. Bullet at lower left: Ac[0] = 0. Left Panel: Ac such 
that Jy = Jz = 0, increasing Jx from to J. Center Panel: Ac 
such that Jx ~ J, Jz ~ 0, increasing Jy from to J. Right 
Panel: Ac such that Jx = Jy ~ J, increasing Jz from to 
J. Connected diamonds: Solution of the full weak/strong 
criterion (^^ forj/hcu — 10. Dashed lines: The estimated 
AsT as given by (|l|) for J/Huj = 1, 2, 5, 10, 100, 1000, from 
top to bottom. 

The abscissa in Figure ^ is essentially the quantity J^/J 
as ranges from to 3 J according to the particular 
scheme chosen for sequentially "turning on" higher di- 
mensions. It is clear that while the casual criterion A '--^ 1 
constitutes a fair order-of-magnitude characterization of 
the occurrence of self-trapping in bulk materials in the 
adiabatic limit, there is considerable qualitative structure 
missed. 

The results shown in Figure ^, though quite general in 
character, depend on the particular manner in which the 
parameters {Jx, Jy, Jz} are varied relative to each other; 
in particular, we note that the manner in which these are 
varied in Figure ^ does not isolate the anisotropy. We can 
obtain a more global view of self-trapping in higher di- 
mensions while simultaneously isolating the anisotropy 
dependence by considering not the composite parameter 
Ac, but the more elementary coupling parameter gc con- 
tained within it according to ([l^). That is. 



9c 



(22) 



This critical value of the coupling constant in the adia- 
batic limit depends on both the intensity of tunneling | J| 
and on the anisotropy. The dependence on the anisotropy 
can be isolated, however, in the normalized quantity 



5c [1] 



1/2 



(23) 
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in which gd^] = '\J\J/hLLj\ represents the critical cou- 
phng parameter in the one-dimensional case subject to 
the condition that the ID tunneling parameter is fixed 
at the value \J\ appropriate to the Z?-dimensional case. 
While dependent on each of Jj;, Jy, and Jz, the ratio 
( p3| ) is independent of | J\ and depends only on the angu- 
lar variables in a spherical polar coordinate representa- 
tion of the {Jx, Jy, Jz} system. Eq. thus describes a 
surface having the interpretation that the radial distance 
from the origin is the factor by which the critical coupling 
constant gc in D dimensions exceeds the critical value in 
one-dimension (5c[l]) having the same intensity of tun- 
neling. This surface is plotted in Figure together with 
an isotropic (spherical) reference surface, and a surface 
corresponding to the condition Ac — const. This latter 
surface shows that the oft cited condition Ac ~ 1 contains 
implicit anisotropy; however, it is evident that the real 
anisotropy of self-trapping is even greater than might be 
inferred from this common rule of thumb. 
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FIG. 4. Dependence of self-trapping on anisotropy. The 
Cartesian axes are the positive half-axes of the {Jx, Jy, Jz} 
system. The innermost surface is an octant of the unit sphere 
exhibited as a reference surface reflecting isotropic depen- 
dence on I J| only. The outermost surface is defined such that 
the radial distance from the origin represents the ratio (?c/ffc[l] 
as in (p^). The interleaved surface is defined such that the 
radial distance from the origin reflects the anisotropy implicit 
in the condition Ac = const, at fixed \J\. 

The presentation of Ac shown in Figure |^ corresponds 
to a particular transit of the g^ surface seen in Figure ^: 
The OD case can be considered to occupy the origin, and 
the ID cases correspond to the three corners of the dis- 
played surface. The "turning on" of the second dimen- 
sion according to the scheme of Figure ^ corresponds to 
movement along the edge of the displayed surface to the 
midpoint of that edge corresponding to the 2D isotropic 



case. The subsequent "turning on" of the third dimen- 
sion corresponds to movement perpendicular from this 
edge along a straight line (geodesic) to the center of the 
surface corresponding to the 3D isotropic case. This com- 
parison shows in particular (as may be proven analyti- 
cally): 1) that both the jump discontinuity between OD 
and higher dimensions and the square-root singularity 
between ID and higher dimensions are generic features, 
not dependent on the manner or sequence with which 
transverse dimensions are "turned on", and 2) that the 
less-singular feature seen in Figure ^ at the transition 
from 2D to 3D is not generic, but appears only because 
dimensions in Figure ^ were turned on sequentially. 

Thus, for a given | J|, we can distinguish three regimes 
based on sensitivity to anisotropy: 

g > gl; there are no large polaron states at any degree 
of anisotropy. 

g < 5c [1]; there are no small polaron states for any 
degree of anisotropy. 

5c[l] < g < gl', large polaron states exist for suf- 
ficiently isotropic tunneling, small polaron states exist 
for sufficiently anisotropic tunneling, and these regimes 
are separated by a self-trapping transition as a function 
of anisotropy at fixed |J| and g. This effect of self- 
trapping as a function of anisotropy alone can be un- 
derstood in terms of the size, shape, and content of the 
phonon cloud. As discussed in Appendix A, the more 
isotropic and higher-dimensional polaron scenarios are 
characterized by phonon clouds that are spread over the 
largest volumes of space and contain the fewest num- 
bers of phonons. With increasing anisotropy, the po- 
laron cloud grows more compressed, occupying smaller 
volumes of space, but being occupied by larger numbers 
of phonons. If this anisotropy-driven compression can 
proceed sufficiently far, the number of phonons in the 
phonon cloud can be driven sufficiently high that self- 
trapping can occur. 

IV. SELF-TRAPPING AWAY FROM THE 
ADIABATIC STRONG-COUPLING LIMIT 

At general parameter values away from extreme limits, 
accurate estimations of the location of the polaron self- 
trapping line are scarce. Until rather recently, estimates 
even in one dimension were largely casual rules of thumb. 
As noted above, a frequently-encountered characteriza- 
tion holds that self-trapping occurs when A ~ 1; this con- 
dition is often supplemented by the condition g > I ac- 
knowledging that the strong-coupling theory from which 
the A condition arises is not expected to hold to arbitrar- 
ily weak coupling. 

We can improve on the common rule of thumb by iden- 
tifying the self-trapping transition not with a single fixed 
value of A (e.g., unity) but with the critical value ob- 
taining in the adiabatic limit for the particular dimen- 
sion and J appropriate to each unique circumstance (i.e., 
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A '-^ Ac). In so doing, we capture all the structure evi- 
dent in Ac (Figures ^ and ^ and make the preliminary 
assumption that the scaling relationships that character- 
ize the adiabatic limit hold to a meaningful degree at 
moderate parameter values; i.e, we may consider extrap- 
olation of critical scaling relationships to finite parameter 
values. The implications of such an assumption for the 
elementary coupling parameter gc are shown in Figure |^. 
The shifting of these estimated self-trapping lines with 
anisotropy is a direct reflection of the anisotropy of Ac- 
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J/hm 



(24) 

This refinement yields estimated self-trapping lines as 
illustrated by the truncated curves in Figure ||; these 
curves are truncated (arbitrarily at J/huj = 2) because 
intersections of WCPT and SCPT begin to disappear at 
lower values of J /huj, as can be seen in the ID panel of 
Figure [I[ 
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FIG. 5. Dependence of the self-trapping line on anisotropy 
under the assumption Ast = Ac- The bold solid lines corre- 
spond to the isotropic cases in (from bottom to top) ID, 2D, 
and 3D. The textured lines between the ID and 2D isotropic 
cases correspond to (from bottom to top) Jx = J , Jz = 0, 
with Jy = 0.2J, 0.4J, 0.6J, and 0.8J. The textured lines 
between the 2D and 3D isotropic cases correspond to (from 
bottom to top) Jx — Jy = J, with = 0.2 J, 0.4J, 0.6J, and 
0.8J. 

The shifting of these estimated self-trapping lines with 
anisotropy is a direct reflection of the anisotropy of Ac; it 
is this qualitative character of the mutual relationships 
among self-trapping curves of differing anisotropies that 
we expect to be largely preserved as necessary correc- 
tions are made. The need for further correction is ev- 
ident, for example, in that the ID example in Figure |^ 
differs substantially in absolute terms from the ID empir- 
ical curve ( |l5| ) although the two are qualitatively quite 
similar. Moreover, all of the curves displayed in Fig- 
ure]^ violate the ancillary condition g > 1 at small J /huj, 
reflecting the expected eventual failure of extrapolation 
from the adiabatic strong-coupling regime. 

We should be able to improve on this estimate by us- 
ing a more complete weak/strong condition ( |2^ ) employ- 
ing the complete results of both perturbation theories 
through second order as given in (^ - (p^, thus ob- 
jectively capturing non-adiabatic corrections implicit in 
those terms that do not contribute in the adiabatic limit. 
Thus we consider the condition 



FIG. 6. Truncated curves show the dependence of 
the self-trapping line on anisotropy using the complete 
weak/strong condition (^^ as a self-trapping criterion; these 
curves are clipped below J/hu = 2 because crossings of weak- 
and strong-coupling results begin to disappear (see the ID 
case of Figure uh- The untruncated solid curves are the em- 
pirical gsT of p5| ) in the isotropic cases of each dimension- 
ality. The bold solid lines correspond to the isotropic cases 
in (from bottom to top) ID, 2D, and 3D. The textured lines 
between the ID and 2D isotropic cases correspond to (from 
bottom to top) Jx ^ J, Jy = 0.2J, 0.4J, 0.6J, and 0.8J, and 
Jz — 0. The textured lines between the 2D and 3D isotropic 
cases correspond to (from bottom to top) Jx = Jy = J, and 
Jz = 0.2J, 0.4J, 0.6J, and 0.8J. 

The effects of including non-adiabatic corrections de- 
pend on dimensionality, anisotropy, and "distance" from 
the adiabatic limit: 1) The self-trapping curves describ- 
ing ID and quasi-lD cases shift strongly to stronger cou- 
pling values, suggesting a corrective shift of order unity 
at essentially all J. 2) The self-trapping curves describing 
2D and 3D cases shift only weakly at moderate adiabatic- 
ity and more weakly with increasing adiabaticity. 3) Ex- 
cept for strong corrections in the quasi-lD regime, the 
qualitative character of the dependence of self-trapping 
on anisotropy is little affected by non-adiabatic correc- 
tions. 4) At low adiabaticity, all self-trapping curves shift 
to stronger coupling values in a manner and to a degree 
consistent with a condition > 1 at J = rather than 
the condition g > suggested by adiabatic scaling. 
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Gathering all the implications of the above together, 
we are led to extend our ID empirical curve (^5|) describ- 
ing the one-dimensional self-trapping to the general case 
describing any dimension and any degree of anisotropy. 
To do this we combine: a) the empirical curve (?st[1] that 
effectually characterizes the one-dimensional case, b) the 
adiabatic critical curve gc that effectually characterizes 
the higher-dimensional, higher-adiabaticity regime, and 
c) the adiabatic critical parameter Ac that compactly de- 
scribes the qualitatively distinct characteristics of the low 
and high dimensionalities. 

From such considerations we are led to a family of em- 
pirical curves 
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(1 + j/M'^=™"'^^-'-^/'-'^'^+5c 



(25) 



in which all quantities have been previously defined. This 
family of curves is not derived from any theory, and, 
apart from the ID case, is not backed by a large body 
of independent high-quality data since such data is quite 
sparse at the present time. What high-quality data does 
exist at the present time is quantitatively consistent with 
this family of curves in the same fashion that an abun- 
dance of high-quality ID data has been found consistent 
with SfsT[l] (see Figures and 0). In keeping with the 
discussion of Section |l| we have not attempted to regu- 
larize square root dependences that arise naturally in the 
adiabatic limit, but which are most likely softened with 
decreasing adiabaticity. The utility of ( |25| ) lies in com- 
pactly and simply describing the apparent and mutually 
consistent trends in a large volume of results of inde- 
pendent methods and arguments, providing meaningful 
estimates for the location of the self-trapping transition 
in any dimension for any degree of anisotropy or adia- 
baticity. This estimated gsT is compared with quantum 
Monte Carlo data for the ground state energy in Figure |l| 
and effective mass in Figure |^. 

In Figure |^, we have included several curves (dashed 
lines) corresponding to \st = d^srl'^J^^ using ( p5| ) for 
J/huj = 1, 2, 5, 10, 100, 1000. These curves indicate 
how we expect the physically meaningful self-trapping 
line at finite parameters as estimated by (|2|) to be re- 
lated to the results of the adiabatic limit. Figure ^ shows 
that the higher-dimensional, more isotropic cases con- 
verge toward their adiabatic limits more rapidly than do 
lower-dimensional, more anisotropic cases. This conver- 
gence in one dimension is particularly poor, with signif- 
icant deviations from the adiabatic limit persisting for 
J/huj > 1000, by which point the higher-dimensional 
cases have converged beyond plotting precision. 

Viewed collectively, the dashed curves of Figure || also 
show that the composite parameter A does not provide a 
very natural or even qualitatively self-consistent charac- 
terization of the self-trapping transition over the whole of 
the adiabatic regime {J/fiui > 1/4). In the far adiabatic 
regime, where we may take Ac to fairly characterize the 
location of the self-trapping transition (solid curve in Fig- 
ure H), one may be led to conclude that large polarons 



are relatively more stable in higher dimensions and at 
weaker anisotropics since the occurrence of self-trapping 
is found to shift to larger values of A in these regimes. 
On the other hand, at more moderate degrees of adia- 
baticity (e.g., J/huj — 1,2,5 in Figure ^ one is led by 
the same reasoning to conclude that large polarons are 
relatively less stable in higher dimensions and at weaker 
anisotropics since the occurrence of self-trapping is found 
to shift to lower values of A in these regimes. In partic- 
ular, one of the most actively-investigated cases in con- 
temporary studies is the "typical" scenario with J/huj of 
order unity; Figure || shows that in terms of A, the self- 
trapping trends in this case are quite distinct from those 
found in the adiabatic limit, certainly complicating the 
interpretation of results. 

From Figure ^, on the other hand, based on the more 
elementary coupling parameter g appearing directly in 
the Hamiltonian, one is led to conclude that large po- 
larons are everywhere relatively more stable in higher 
dimensions and at weaker anisotropics since the self- 
trapping line shifts to larger values of g as these trends are 
followed regardless of the degree of adiabaticity. These 
trends in g are qualitatively similar and uniform for all 
degrees of anisotropy and adiabaticity, whereas the same 
trends in A vary strongly with regime. For the same rea- 
sons that A is a convenient parameter with which to char- 
acterize polarons in the far adiabatic regime, it proves to 
be an inconvenient parameter in the broader context of 
the problem away from the adiabatic limit. 



V. CORRELATION FUNCTION AND POLARON 
RADIUS 

In view of the local nature of the electron-phonon cou- 
pling in the Holstein model, the spatial extent of the 
polaron can be characterized quite directly through an 
analysis of electron-phonon correlations. This can be 
done using a correlation function that has been long and 
widely used to characterize polaron size in D dimensions 
©Hill: 



= (C-f ) = ^ ^(al-a,,(6l.+, + bn+,)) , (26) 

n 

normalized such that J^r^l^^ = 1- This function can 
be viewed as measuring the shape of the polaron lattice 
distortion around the instantaneous position of the elec- 
tron. 

Using Rayleigh-Schrodinger perturbation theory in the 
weak-coup ling re gime as in the preceeding sections one 
finds that ||27| , |35|l 

Ct°i ^hcu dt e-^^'T\e-^J^'lr^{2J,t) . (27) 
•^0 z=l 

Note that setting any one Ji to zero or summing ct^' over 
one ri recovers ct^ ^' . This property implies that the 
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effect of "turning on" transverse dimensions is simply to 
spread electron-phonon correlation strength transversely. 

Characterizing this multi-dimensional correlation func- 
tion in terms of a width measure involves a variance ten- 



sor, 



Contrary to much prevailing opinion, these results 
show that in the weak-coupling regime: i) there are no 
significant qualitative or quantitative differences between 
ID, 2D, and 3D polaron radii, ii) the polaron radius in 
2D and 3D is not infinite, and iii) the polaron radius 
does not scale as J / g^TiLo in any dimension as commonly 



(28) expected, but as J jfiuj in every dimension 



where 



/>oo 

<jl = hLj dt e-''^' rfe'^'^^*Ir^{2J,t) (29) 
Jo 

(30) 



2J, 

huj 



h 1 

2m?,w 1? 



in which 



is the free electron effective mass and /, is 



the lattice constant in the i direction. Thus, along each 
of the primitive crystallographic axes, the real-space vari- 
ance is simply proportional to the electron transfer inte- 
gral along that axis, and in a general direction is just 
the appropriate mixture determined by rotation. In ab- 
solute units (unrationalized by the lattice constants) the 
real-space variance is the same as that of the zero-point 
motion of a harmonic oscillator characterized by the lat- 
tice frequency w, but with the lattice mass replaced by 
the free electron mass measured along the appropriate 
direction. 

Utilizing the notion of a polaron half-width defined in 
terms of the correlation variance 



Ri 



(31) 



we can associate with the polaron characteristic ellip- 
soidal volumes V[D] 



2Jx 



1/2 



2uj 



(32) 



1/2 



V[3] 



4 V2w 
47r 



- — (det M, 



-1n1/2 



(33) 



ODD 1 1 1 ^ f ^"^^ '^'^y ^"^^ 



1/2 



TT f h 



6 V2w 



3/2 



[det Mo^) 



IN 1/2 



(34) 



This characteristic volume thus increases with the inten- 
sity of tunneling (V^[I?] oc \J /%uj\^/'^)^ and is largest 
in the isotropic case and decreases with increasing 
anisotropy. In the isotropic case we may regard R = Ri 
as the polaron radius. 



VI. EFFECTIVE MASS 

For the circumstances we address in this paper, the 
reciprocal effective mass tensor is diagonal, with elements 
given by 



{M- 



dKidKi 



^6.~. (35) 



K = 



From this, it is easily shown that the reciprocal effec- 
tive mass in any direction through second order of weak- 
coupling perturbation theory is given by 

D 



rUi. 



/>oo ^ 

l-g^h'^oj^ dt t e-'^"*TTe-2-^.*/o(2J,t), (36) 
•^0 ^=l 



where m*j and are respectively the polaron and free 
electron effective masses in the i direction.. 

Figure ^ shows the dependence of the isotropic polaron 
mass on dimensionality according to WCPT and quan- 
tum Monte Carlo simulation. Although this is a compar- 
ison between isotropic cases of J/Tllu — 1 only, the ex- 
cellent agreement between WCPT and quantum Monte 
Carlo out to ~ \/D suggests that the WCPT mass may 
be similarly accurate for A < 1/2 as defined in (|l^ ) at 
general J. 




0.6 



0.2 



WCPT 

Global-Local 

a— CQMC 

▼ STT 




0.0 



FIG. 7. The reciprocal effective mass ratio m%/m*i for 
the isotropic case. Solid curve: Global-Local mass for D = 1 
Chain-dashed curves: WCPT masses for D = 1, 2, and 3 
Scatter-plot: Quantum Monte Carlo data for D = 1, 2, and 3 
data kindly provided by P. E. Kornilovitch pUp^ . Triangles 
Estimated self-trapping points gsT as discussed in Section [V 
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The weak-coupling result (jS^) shows that although 
anisotropy has definite effects on the value of the effective 
mass, the effect of anisotropy appears only in the value 
of a scalar multipler of the free electron mass; that is, 
although anisotropy of the free electron mass (inequal- 
ities among J^, Jy, and J^) is manifested in real-space 
anisotropics in electron-phonon correlation (i.e., in dis- 
tortions of the shape of the polaron as discussed in the 
previous section), the mass renormalization associated 
with such distortions of polaron shape is isotropic. In- 
terestingly, this implies that increasing Jy or at fixed 
Jx (for example) results in a decrease in m.*^, translating 
into an associated increase in mobility in the x direc- 
tion. This influence of transverse directions on to*^. is 
illustrated in Figure ^. In the center and right panels of 
Figure ^, Jxj'huj is held fixed at unity, yet the effective 
mass in the x direction continues to decrease as tunneling 
into transverse dimensions is turned on. 

These effects can be understood in terms of the trans- 
verse spreading of electron-phonon correlation strength 
as discussed in the last section. As a fixed correlation 
strength is spread over an increasing number of sites 
(characteristic volume of the polaron increases as dis- 
cussed in the previous section), the average lattice defor- 
mation per participating site decreases. Consequently, 
mean square measures of lattice deformation decrease 
and exhibit changes that suggest a diminishing effective- 
ness of electron-phonon interactions in producing typical 
polaronic effects. The polaronic mass enhancement bears 
such a mean square dependence on the lattice deforma- 
tion and like other such measures (e.g., the number of 
phonons in the phonon cloud as discussed in Appendix 
A) , decreases with increasing dimensionality and decreas- 
ing anisotropy. 
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FIG. 8. The effective mass ratio m'^/m^j, according to 
(p6[), showing the effects of changing anisotropy and dimen- 
sionality. Left Panel: ID mass (Jy = Jz ~ 0) at fixed 
g, increasing Jx/ftu from to 1. Center Panel: 2D mass 
{Jx/hu, Jz/hu! = 0) at fixed g, increasing Jy/hu) from to 
1. Right Panel: 3D mass [Jx/Tiuj = Jy/Tiuj = 1) at fixed g, 
increasing Jz/fioJ from to 1. 

The corresponding polaron effective mass resulting 
from strong-coupling perturbation theory through sec- 
ond order is given by 

!4 = ^-9' + e-25V(5')(3J. + J)/noj . (37) 

This result is isotropic at first order simply by virtue of 
being independent of all Ji at that order, but the second- 
order correction is anisotropic because the r.h.s. of ( |37| ) 
bears an explicit, unbalanced sensitivity to the direction 
along which the effective mass component is being mea- 
sured. 

Unfortunately, this strong-coupling result is not very 
helpful; it disagrees substantially with more reliable re- 
sults pO|- p^ , ^ except at small J/Tiuj. We take this as 
an indication that dominating (perhaps non-exponential) 
contributions have yet to be extracted from higher orders 
of SCPT. For such reasons we cannot estimate the loca- 
tion of the self-trapping transition from any crossing of 
( ^ ) and (|37|). Instead, we have included in Figure]^ sev- 
eral symbols to indicate the values of gsr as given by 
); these several values are mutually consistent in lo- 
cating essentially the same feature of the effective mass 
in every dimension, and coincides very well with the ef- 
fective mass feature we have previously identified with 
the self-trapping transition (see Ref. ||3l|]). 

VII. ON DIMENSIONALITY AND 
ADIABATICITY 

As noted in the introduction, the results that have 
long characterized commonly-held expectations for the 
dimensionality dependence of polaron structure arc due 
to behavior ascertainable in the adiabatic approximation 

In 2D and 3D, the minimum energy states in the 
adiabatic approximation are found to be "free" states 
throughout the weak-coupling regime up to a discrete 
coupling threshold beyond which "self-trapped" states 
have the minimum energy. This abrupt transition phe- 
nomenon is what is meant by the term "self-trapping 
transition" in the adiabatic approximation. Accordingly, 
there is no occasion to distinguish large polarons from 
small polarons in 2D and 3D since the "free" states below 
the transition are of infinite radius and distinct from large 
polarons, and the "self-trapped" states above the transi- 
tion are always interpretable as small polarons. This set 
of circumstances in 2D and 3D is reflected in the catch 
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phrase "all polarons are small", since in this view large 
polarons in the adiabatic sense are never characteristic 
of the polaron ground state in bulk materials. 

In ID, on the other hand, "free" states are unsta- 
ble in the adiabatic approximation; instead, finite-radius 
(i.e. "self-trapped") states are found at all finite cou- 
pling strengths, leading to the commonly encountered 
view that there is no self-trapping transition in ID. That 
polaron states in ID might be distinguishable as large or 
small is inconsequential in this view, as is the notion of 
a resolvable transition between distinct large and small 
polaron structures. 

The results of this paper differ strongly from the con- 
ventional adiabatic picture in multiple respects: 

i) The quasiparticles implicit in the weak-coupling 
states of every dimension are not weakly-scattered "free" 
electrons, but dressed electrons having finite radii gener- 
ally greater than a lattice constant. 

ii) Although these weak-coupling quasiparticles can be 
sensibly characterized as large polarons, in no dimension 
do these weak-coupling states coincide with the large po- 
laron states familiar from the adiabatic approximation in 
ID. 

iii) The finite radii characterizing the weak-coupling 
quasiparticles in every dimension saturate to finite values 
with vanishing electron-phonon coupling, unlike the large 
polaron radii in the adiabatic approximation that in ID 
diverge with vanishing coupling and in 2D and 3D are 
infinite already at finite coupling. 

iv) The self-trapping transition exists in every dimen- 
sion, including ID. 

v) The self-trapping transition is associated with the 
change from large polaron structure to small polaron 
structure in every dimension, including 2D and 3D, and 
not with a change from infinite to finite radii. 

vi) Dependences of polaron properties on parameters 
are smooth through the self-trapping transition in every 
dimension, unlike the abrupt changes often found in the 
adiabatic approximation in 2D and 3D. 

Our results are quantitatively supported by indepen- 
dent high-quality methods (including variational meth- 
ods [|^]3l|,||,|^,|3, cluster diagonalization 
density matrix rcnormalization group 1 20 1 , and quantum 
Monte Carlo (l|-|l|j2l|,|2|,||l ) . Moreover, elaborations 
of adiabatic theory incorporating non- adiabatic correc- 
tions support our overall conclusion that the 
adiabatic approximation as it is widely regarded fails to 
embrace non-adiabatic characteristics that are essential 
to the proper description of polaron states in the weak 
coupling regime, and therefore fails as well to properly 
describe the self-trapping transition itself Q]. 

With so many results at variance with the adiabatic 
approximation, it is well to ask in what respects, if any, 
are our results consistent with the adiabatic approxima- 
tion and whether some sense can be made of the pervasive 
discrepancies. Indeed, several consistencies can be found 
that are illuminating. 



We first note that the dependence of Ac on dimen- 
sionality and anisotropy exhibits a generic square-root 
singularity at the boundary between ID and any higher 
dimensional case, while at the boundary between higher- 
dimensional cases this dependence is generically smooth. 
For essentially the same underlying reasons, Ac is con- 
stant throughout ID, but varies with detail of tunnel- 
ing in higher dimensions. These distinctions are at least 
suggestive of the sharp contrasts between ID and higher 
dimensional cases in the adiabatic approximation. 

Secondly, we note that the weak-coupling polaron ra- 
dius R as here derived diverges in any dimension in the 
adiabatic limit. Further considering the WCPT validity 
test in Appendix A, there is reason to speculate that this 
weak-coupling radius might continue to be a reasonably 
valid construct in 2D and 3D up to the vicinity of the 
self-trapping transition. Such a possibility might be con- 
sistent with the finding of strictly infinite-radius states 
on the weak-coupling side of the transition in 2D and 
3D in the adiabatic approximation, while the possible 
breakdown of the weak-coupling radius construct below 
the transition in ID might be consistent with existence 
of finite-width states in the ID adiabatic approximation. 



VIII. CONCLUSION 

In this paper we have analyzed the dependence of nu- 
merous polaron properties on the effective real-space di- 
mensionality and anisotropy as determined by the elec- 
tronic tunnelling matrix elements; these properties in- 
clude the polaron ground state energy, polaron shape, 
size, and volume, the number of phonons in the phonon 
cloud and the polaron effective mass. In pursuing 
these analyses we have made extensive use of weak- 
and strong-coupling perturbation theories supported by 
selected comparisons with non-perturbative methods. 
Through the use of a scaling argument combining weak- 
and strong-coupling perturbation theory in the adiabatic 
strong-coupling regime, we have been able to infer the 
probable location of the self-trapping critical point in the 
adiabatic limit in any dimension and for any degree of 
anisotropy, and by combining information from multiple 
sources we have been able to extend this estimate from 
the adiabatic limit to finite adiabaticity. 

Central among our findings is the over-arching qualita- 
tive conclusion that polarons in any dimension and any 
degree of anisotropy are similar in most respects. In par- 
ticular, polarons on the weak-coupling side of the self- 
trapping transition share a structure that is essentially 
identical in every dimension. This weak-coupling struc- 
ture is consistent with the notion of the weak-coupling 
polaron as a finite-radius quasiparticle, but is incon- 
sistent both with the notion of a weakly-scattered free 
electron (adiabatic approximation in 2D and 3D) and 
with the historical notion of the large polaron (adiabatic 
approximation in ID). The strong-coupling structure is 
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consistent with traditional notions of small polarons, in- 
cluding strong-coupling perturbation theory and the adi- 
abatic approximation. 

Since the essential character of the weak-coupling 
states and strong-coupling states is only inessentially af- 
fected by dimensionality and anisotropy, the notion of the 
self-trapping transition separating the weak- and strong- 
coupling states is similarly not altered in any essential 
way by changes in dimensionality or anisotropy. Neces- 
sarily, one is led to view self-trapping as the more-or-less 
rapid transition, occurring in every dimension, between 
characteristic weak- and strong-coupling states, both of 
which are characterized by finite radii. 

If we may transcend the jargon that historically has 
had a tendency to polarize the conventional wisdom, it 
is fairly concluded that not all polarons are small, even 
in bulk materials, and that in every dimension and for 
every degree of anisotropy the self-trapping transition is 
a smooth, albeit rapid crossover between large and small 
polaron character. 
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APPENDIX A: BREAKDOWN OF 
WEAK-COUPLING PERTURBATION THEORY 

The weak-coupling perturbation theory considered in 
this paper is based on an expansion in states contain- 
ing limited numbers of phonon quanta. The zeroth or- 
der properties are based upon states containing zero 
phonons, and second order properties upon states con- 
taining one phonon. The first neglected order of WCPT 
is the fourth order, built upon states containing no more 
than two phonons. A test of internal consistency of 
WCPT at particular parameters, therefore, is to com- 
pute the expected number of phonons to the retained 
order of perturbation theory, and compare this number 
to the maximum number of phonons present at that or- 
der. For second order WCPT as used in this paper, this 
number of phonons should be less than unity. 

The required computation is contained in 
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relative to huj, one finds that 
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(A3) 
(A4) 



These expressions can be consolidated into the single ap- 
proximate relation 



nph oc g 



MP] 

V[D] 



(A5) 



where il[D] is the primitive cell volume and V[D] the 
characteristic volume of the polaron in D dimensions. 
The dimension-dependent constant of proportionality is 
near 1/2 in all cases. This simple relation, here proven 
only for the adiabatic weak-coupling regime (broad po- 
larons), demonstrates the very direct but inverse relation 
between the number of phonons in the phonon cloud and 
the volume occupied by it. 

In the isotropic case and in terms of the composite 
parameter A, the condition that expected phonon num- 
bers should be less than unity results in the conditions 
(J > hcu) 



A<2 (-^ 

Tiu! 
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< 
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in ID , 
in 2D , 

in 3D . 



(A6) 
(A7) 

(A8) 



Recalling that the self-trapping transition is expected 
to occur at A of order unity, it would appear that WCPT 
through second order is consistent with the condition 
rip/i < 1 up to the transition in 2D and beyond the tran- 
sition in 3D. It is the ID case that appears to be on the 
weakest footing in the adiabatic strong-coupling regime; 
however, it is the ID case that has been most exhaus- 
tively studied by non-perturbative means and found to 
be widely consistent with second-order WCPT. 
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